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findings are that subsonic aircraft emissions: 1) have 
not be responsible for the observed water vapor trends 
at 40°N; 2) could be a significant source of soot mass 
near 12 km, but not at 20 km. 3) might cause a notice- 
able increase in the background sulfate aerosol surface 
area and number densities (but not mass density) near 
the northern mid-latitude tropopause, and 4) could pro- 
vide a global, annual mean top of the atmosphere ra- 
diative forcing up to -f 0.006 W/m 2 and -0.013 W/m 2 
due to emitted soot and sulfur, respectively. 

Introduction 

Atmospheric effects of aviation are a focus of ongo- 
ing NASA and European Programs [e.g., Stolarski et 
al , 1995; Friedl, 1997; Schumann et al, 1997; Brasseur 
et al, 1998] and the forthcoming IPCC Special Re- 
port on Aviation and the Global Atmosphere. Until re- 
cently, most research concentrated on the ozone layer’s 
response to emissions of NO x and H 2 O by the current 
subsonic and projected supersonic aircraft fleets. How- 
ever, in situ measurements of the Concorde exhaust [Fa- 
hey et al, 1995] showed a surprisingly large amount of 
very small sulfate particles, presumably due to rapid ox- 
idation of exhaust SO 2 * Subsequent model calculations 
[ Weisenstein et al, 1996] showed these small particles, 
despite their small total mass, could significantly en- 
hance the stratospheric surface area density (SAD) and 
alter the ozone balance in the lower stratosphere. Like- 
wise, the recent SUCCESS campaign studied the poten- 
tial for aircraft to produce persistent contrails and pos- 
sibly to induce cirrus formation [ GRL , nos. 8-10, 1998]. 

Thus aircraft emission of particles (sulfate, soot, met- 
als, etc.), or of particle precursors (SO 2 , H 2 O), have 
become a significant source of uncertainty in current 
environmental assessment studies. 

Evaluating the environmental impact of the subsonic 
fleet is a very difficult issue due to the following prob- 
lems. First, emissions occur near the tropopause, 
which is a region of large gradients in most chemical 
species and is the interface between two very differ- 
ent meteorological regimes; stratosphere/troposphere 
exchange is poorly understood: and measurements near 
the tropopause are as difficult to interpret as to model. 

Second, the transition in aerosol chemistry across the 
tropopause is complex as it changes from liquid sulfu- 
ric acid and water to a more complicated liquid-solid 
mix in the troposphere that includes ammonia, organ- 
ics, and crustal materials [GRL, no. 9, 1998]. Third, 
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aerosol distributions in the UT/LS are inhomogeneous, 
especially near convective regions or flight corridors, re- 
quiring computationally intensive 3-D models to sim- 
ulate their patterns rather than the longitudinally av- 
eraged 2-D models traditionally used for stratospheric 
ozone assessments. 

This study uses a suite of global 3-D and 2-D models 
to predict the atmospheric accumulation and dispersion 
of exhaust products from the 1992 subsonic fleet. We 
intercompare the models to derive uncertainty ranges 
since no measurements exist that uniquely allow us to 
measure the accumulation of subsonic exhaust products 
in the UT/LS region. Upper limits for the buildup of 
aviation-induced aerosols as well as stratospheric H 2 0 
are derived. We consider whether subsonic emissions 
could explain recent measurements of stratospheric soot 
or water vapor. The impacts on heterogeneous chem- 
istry through changes in surface area density are also 
evaluated, as well as the changes in radiative forcing of 
climate through aerosol-radiation coupling. 

Aviation Fuel Tracer Simulation 

Four 2-D and seven 3-D global model choose to par- 
ticipate in this study (see Table 1). We assume that 
these eleven models can provide uncertainty bounds for 
the model calculations. The 3-D models have indepen- 
dently derived meteorologies. Some focus on the tro- 
posphere and lower stratosphere (e.g. ECHAM3, UIO, 
TM3), while others cover the whole stratosphere as well. 
The 2-D models have independently formulated zonal- 
mean residual transport and diffusion and focus on sim- 
ulating stratospheric ozone. For this fuel tracer sim- 
ulation, each model provided identical diagnostics on 
a standard grid (o° latitude by 2 km altitude) of the 
steady-state, zonally and annually averaged distribu- 
tion of aviation fuel emitted according to a 3-D pattern 
of flight routes designated in the 1992 NASA scenario 
[Baughcum et al , 1996] (Figure 1, bottom right, 139.4 
Mt/vr). The fuel is transported in the models as a pas- 
sive tracer. The only tracer sink was globally uniform 
removal everywhere below 400 hPa (~7 km) with an 
e-folding time of 5 days, simulating tropospheric wash- 
out. 

The mass distribution of a specific exhaust product, 
X, that is reasonably simulated by this removal process 
can be obtained by multiplying the tracer mass mixing 
ratio by the emission index (El) of X, in kg(X)/kg(fuelj 
for the purposes of this study. The following Els are 
used below: EI(H*20) = 1.23, El(soot) =4xl0“ 5 , and 
EI(S) = 4xl0 -4 kg/kg. Volume mixing ratios for gases 
can De derived from scaling the mass mixing ratio by 
/lair! ^tracer (here fj, is the molar weight). 

Results and Discussion 

Figure 1 and Table 1 summarize the principal results 
of the fuel tracer simulation. The maximum in tracer 
mixing ratio occurs for all but the TM3 model between 
50°N and 65°N near 10-12 km, consistent with the 
source distribution. The maximum tracer mixing ratio 
ranges from 12.6 ng/g (ECHAM3) to 122 ng/g (GSFC- 
2D), reflecting in part the model resolution and in part 
the uncertainty in modeling the dispersion of aircraft 
exhaust (i.e. the GSFC-2D and LLNL models have the 
smallest K :z and K yy values there, while the ECHAM3 


model has the strongest vertical transport). To smooth 
the effects of model resolution we calculate the tracer 
amount in the 8-16 km altitude region between 30°N 
and 90°N (shown by the thick dashed line in Figure 1). 
This region contains from 34% (AER, ECHAM3) to 
61% (GSFC-2D) of the total accumulated tracer. The 
absolute amount of tracer mass in this volume ranges 
from 2.9 Tg (ECHAM3) to 14.5 Tg (GSFC-2D). The 
compact latitudinal containment of the tracer poleward 
of 30°N is an encouraging similarity among all models. 

The amount of tracer transported upward towards 
the stratospheric ozone maximum is an important di- 
agnostics in the model results. The amount of fuel 
tracer above 12 km ranges from 14-16% (UCI/GISS 
and GSFC-2D) to 44-45% (UMICH and AER) of each 
model’s global tracer amount. The models also show 
very different tracer mass mixing ratios above the lower- 
most stratosphere, reflecting uncertainties in modeling 
upward motions and diffusion out of the mid-latitude 
tropopause region. For example, at 26 km and 60°N 
the results in Figure 1 vary from 0.9 (UCI/GISS) to 
11.2 ng/g (UMICH). The global residence time of the 
fuel tracer, defined as the ratio of the steady-state tracer 
mass to the tracer source, varies from 21-22 days (TM3 
and ECHAM3) to 62-65 days (GSFC-2D and LLNL). 

H 2 O. The results of this simulation are used to es- 
timate upper bounds for the contribution of subsonic 
aviation to the abundances and trends of H 2 O by scal- 
ing the tracer distribution by the EI(H 2 0). A caveat 
is that the simple uniform removal timescale of 5 days 
below 400 hPa may not represent the removal of ex- 
cess water vapor by the hydrologic cycle. Among the 
models, the LLNL model gives the largest tracer accu- 
mulation at 40°N in the lower stratosphere, with values 
gradually decreasing from 34 ng/g (or 55 ppbv) at 10 
km to 8 ng/g (or 12 ppbv) at 24 km. Assuming 5%/yr 
growth in fuel consumption and EI(H 2 0)— 1.23, one de- 
rives a trend of the aircraft-produced H 2 0 ranging from 
3.4 ppbv/yr at 10 km to 0.8 ppbv/yr at 24 km. For am- 
bient values of water vapor at 40° N of 59 ppmv at 10-12 
km and 4.2 ppmv at 22-24 km [Oltmans and Hofmann, 
1995], the growth in air traffic could produce a trend 
of +0.006 %/yr at 10 km and +0.018 %/yr at 24 km. 
These trends are a factor of more than 20 less than those 
reported by Oltmans and Hofmann [1995]. 

Soot. Calculations for soot with El(soot) = 4x10-5 
kg/kg show the strongest local perturbation of 1.8 
ng/m 3 in the GSFC-2D tracer field at 10 km and 55°N. 
Again, caveats must be made regarding the realism of 
atmospheric scavenging of soot by the parameterized 
loss below 400 hPa. This zonal mean soot mass den- 
sity is comparable to measured values of up to 2 ng/m 3 
at 10-12 km [Pueschel et al, 1997] and recent model 
calculations [Rahmes et al., 1998], indicating that sub- 
sonic aircraft can be a noticeable source of soot mass 
near the northern mid-latitude tropopause. This con- 
clusion is uncertain because we took the maximum 
predicted concentration and existing measurements are 
sparse [Pueschel et al, 1997]. Some other 3-D model 
calculations [Cooke and Wilson , 1996; Liousse et al, 
19961 indicate the importance of surface soot sources 
at 10-12 km altitude. Aircraft, however, are unable to 
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explain the observed soot at 20 km [Pueschel et al, 
1997). The fuel tracer simulation does provide a basis 
to estimate of the perturbation of soot SAD, which is 
a necessary parameter for heterogeneous chemistry cal- 
culations. For example, from the mass balance and as- 
suming the radius of particles r=10 nm, density p 300i - 
2 g/cm 3 , EI(soot)=4x 10 -5 kg/kg, p air = 0.38 kg/m 3 , 

one gets SAD soot = 0.28 pm 2 /cm 3 for the maximum of 
the GSFC 2-D output. 

Sulfate aerosol. The LLXL model results show' the 
largest amount of sulfate aerosol mass in the strato- 
sphere (16.6 ktons, or only 2.1% of its background value 
of 0.8 ktons [G.K. Yue, personal communication]) as- 
suming EI(S)=4xlO~ 4 kg/kg and 100% conversion of 
sulfur emissions to sulfate aerosol. As a result, the trend 
in the stratospheric sulfate aerosol mass due to the sub- 
sonic fleet is 0.1 %/yr if the fuel consumption grow'th 
rate is 5%/yr and natural sources are constant. This 
trend is much smaller than the 5%/yr trend found by 
Hofmann [1991] and contradicts his inference that the 
observed trend was due to the sulfur in aircraft fuel. 
However, our conclusions are consistent with those in 
Bekki and Pyle [19921. 

The contribution of subsonic aircraft to the sulfate 
aerosol mass at 10 km and 55°N is 110 ng/m 3 (GSFC- 

2D) for EI(S) of 4xl0“ 4 kg/kg and 100% conversion. 
This value is several times smaller than the background 
value [Yue et al . , 1994]. Thus, even at the point of max- 
imum accumulation, subsonic aircraft could not perturb 
the sulfate aerosol mass strongly. 

The rapid formation of new nm-sized sulfate particles 
in aircraft plumes is estimated to vary betw*een 0.4% 
[Curtius et al, 1998] and >12% [Fahey et al, 1995] of 
available sulfur emissions. The new' sulfate particle for- 
mation in the plume determines the perturbations of 
SAD and hence the heterogeneous chemistry [Weisen- 
stein et al, 1996]. The remaining S0 2 gas will oxidize 
more slowly and condense onto existing sulfate parti- 
cles, thus contributing to the visible scattering and, 
hence, climatic impacts. Aircraft perturbations to the 
sulfate aerosol SAD at altitude z and latitude 0 are es- 
timated as: 


SAD su |(z, 0) = l t .!;( Z ^)Pair^E I(S)^H2S04 
rp su |WtH2S04^S 


(i) 


Here, tr(z, 0) is the tracer distribution, f is the con- 
version fraction of sulfur to sulfate, and wtH2SO\ is the 
weight fraction of H 2 S0 4 in a particle. Since the un- 
certainties of the aircraft sulfur emission conversion to 
new sulfate particles are large, we assume £=5% for 
demonstration purposes only. For example, assuming 
r= 10 nm, p air = 3.8xl0" 4 g/cm 3 , p 3ul ~ 1.5 g/cm 3 , 
EI(S) =0.0004, and wt// 2 so4 := 0.5 (or 50%), equation (1) 
gives 1.13 /im 2 /cm 3 for the GSFC-2D maximum of 122 
ng/g. This value is comparable to background sulfate 
SAD and consistent with the results of Kjellstrom et 
al [19981; however, it becomes much smaller than back- 
ground levels at 16-20 km where heterogeneous chem- 
istry plays an important role in mid-latitude ozone loss. 

Direct Radiative Forcing of Climate. The atmo- 
spheric accumulation of the fuel tracer is non-uniform 
with maximum zonal annual mean concentrations of 



4.1-23 /ig/cm 2 located along the flight corridors near 
50°-60°N. The global mean columns range from 1.6 to 
5.0 jxg/cm 2 with a median value of 2.3 jxg/cm 2 . The 
potential climate change caused by aircraft is assessed 
here by the global, annual mean instantaneous top of 
atmosphere radiative forcing (TOARF). 

The global mean column burden of soot ranges from 
0.07-0.20 ng(soot)/cm 2 . The best current estimate 
of the sensitivity of TOARF to soot column burden, 
~3000 W/g(soot), includes enhanced effects of absorb- 
ing aerosols in the upper troposphere above clouds 
[Haywood and Ramaswamy, 1998]. If aircraft-produced 
and background soot have similar properties, an upper 
limit of TOARF (soot) is +0.006 W/m 2 . 

Assuming an EI(S) of 4xl0 -4 and 100% conversion 
of sulfur to sulfate, the global mean S0 4 column bur- 
den is 2. 0-6.0 ng(S0 4 )/cm 2 . Assuming a sensitivity 
of -215 W/g(S0 4 ) for relatively dry sulfate aerosols in 
the UT/LS [Haywood and Ramaswamy , 1998], the up- 
per limit of the TOARF(sulfate) is -0.013 W/m 2 . Our 
TOARFs provide upper limits, since coagulation and 
sedimentation are ignored in our tracer simulation and 
100% conversion of SO 2 sulfur is assumed. Also, models 
with aerosol microphysics [Chuang et al, 1997] suggest 
~50% conversion. These TOARFs are approximately 
two orders of magnitude smaller than that due to cu- 
mulative anthropogenic CO 2 emissions from fossil fuel 
and are of similar magnitude to TOARF due to subsonic 
aircraft CO 2 emissions. 

Concluding remarks. Our results from a tracer 
study for the 1992 aircraft fleet provide a useful quan- 
titative assessment of aircraft effects for a wide range 
of potential atmospheric perturbations. The indepen- 
dent models participating show a modest range for the 
global mean accumulation of aircraft exhaust, ranging 
only over a factor of three despite the large differences in 
modeling the tropopause region. Our results are most 
useful as upper bounds to the accumulation of aircraft 
exhaust products because our simulation ignores pho- 
tochemical and microphysical processes. Despite these 
limitations, we are able to place the aircraft environ- 
mental effects in perspective and rule out some hypothe- 
ses attributing observed changes in sulfate aerosol mass 
and H 2 O to aircraft. More accurate assessments of the 
atmospheric effects of aircraft will require 3-D atmo- 
spheric models with full photochemistry and/or micro- 
physics. 
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Figure 1 . Zonally and annually averaged distribution 
of the fuel tracer in ng(tracer)/g(air) from the partici- 
pating models. The fuel source is shown in the bottom 
right panel. The thick dashed line shows the region 
between 8-16 km and 30°N-90 o N, the thin dashed line 
depicts the 12 km altitude. 

Figure 1. Zonally and annually averaged distribution of the fuel tracer in ng(tracer)/g(air) from the participating 
models. The fuel source is shown in the bottom right panel. The thick dashed line shows the region between 8-16 
km and 30°N-90°N, the thin dashed line depicts the 12 km altitude. 
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Table 1. Results from the 1992 aviation fuel tracer simulation 


Model 1 

Resolution: 

Maxi- 

Lat. 

Global 

Tracer in 

% of 

Max. 3 

Global 

Global -1 

Model 5 


lat xalt 

mum. 

of 

resid. 

8- 16km 

tracer 

tracer 

soot 

S0 4 

descri- 


or 

lat x Ion xalt 

ng/g 

max. 

°N 

time, 

days 

30N-90X, 

% 

above 
12 km 

column 

/ig/cm 2 

column 

ng/cm 2 

column 

ng/cm 2 

ption 


AER 

9.5°xl.2km 

26.7 

55 

38 

34 

45 

6.6 

0.11 

3.5 

W98 

GSFC-2D 

10°x2km 

122 

55 

62 

61 

16 

22.9 

0.20 

5.9 

J96 

LLNL 

5°xl.5km 

72.5 

65 

65 

42 

38 

14.5 

0.20 

6.0 

K94 

UNIVAQ-2D 

10°x2.8km 

36.4 

60 

23 

58 

33 

7.7 

0.08 

2.2 

P93 

ECHAM3 

5.6°x5.6°x2km 

12.6 

50 

22 

34 

31 

4.1 

0.07 

2.0 

SK94 

GSFC-3D 

2°x2.5°x2km 

46.7 

50 

52 

44 

29 

11.7 

0.16 

4.9 

W96 

TM3 

3.75° x5°xl. 5km 

20.1 

80 

21 

45 

40 

4.9 

0.07 

2.0 

Wa97 

UCI/GISS 

4 0 x5°x2.5km 

34.4 

55 

27 

49 

14 

8.2 

0.09 

2.6 

H98 

UIO 

8°xl0°x3km 

28.2 

55 

29 

50 

40 

7.8 

0.09 

2.7 

BI97 

UMICH 

4°x5°x2.75km 

30.4 

65 

45 

37 

44 

9.6 

0.14 

4.2 

P91 

UNIVAQ-3D 

10°x22.5°x2.8km 

38.4 

50 

25 

50 

41 

7.7 

0.08 

2.3 

Pi93 


l The models are denoted as follows: (2-D) Atmospheric and Enironmental Research, Inc (AER), Goddard Space Flight 
Center (GSFC-2D), Lawrence Livermore National Laboratory (LLNL), University of L’Aquila (UNIVAQ-2D) and (3-D) 
DLR application of Hamburg GCM (ECHAM3), GSFC-3D model with cloud mass fluxes, Royal Netherlands Meteorological 
Institute and Utrecht University (TM3), University of Calfornia at Irvine (UCI/GISS), University of Oslo (UIO), University 
of Michigan (UMICH), and UNIVAQ-3D. 

2 Vertical resolution shown corresponds to the levels near the tropopause. 

3 Column amounts are calculated from 0 to 60 km for all models except ECHAM3 and TM3 (up to 32 km) and UIO (up 
to 26 km). 

4 These values assume EI(sulfur)=4x 10~ 4 kg(S)/kg(fuel) with 100% conversion of sulfur to sulfate particles. 

5 W98: Weisenstein et a/. [1998]; J96: Jackman et al. [1996]; K94: Kinnison et ai [1994]; P93: Pitari et al. [1993]; SK94: 
Sausen and Kohler [1994]; W96: Weaver et al. [1996]; Wa97: Wauben et al. [1997]; H98: Hannegan et al. [1998]; BI97: 
Bemtsen and Isaksen [1997]; P91: Penner et al. [1991]; Pi93: Pitari [1993]. 
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